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ABSTRACT 

We report on ASCA energy spectra of Cir X-l taken near its zero phase on 1994 
August 4-5. The ASCA SIS detectors allow a much more detailed study of the iron 
K complex than has been possible before. We find that prior to a sudden upward flux 
transition the dominant iron K feature appears to consist of a large edge from neutral 
or nearly-neutral iron. The depth of the edge corresponds to an absorption column of 
« 1.5 x 10 24 cm~ 2 , while little absorption over that expected from the Galaxy is seen 
at lower X-ray energies. The differential absorption at high and low X-ray energies 
combined with the iron edge energy are strong evidence that partial covering is a 
crucial determinant of the behaviour observed from Cir X-l. The continuum spectral 
variability observed by ASCA can also be understood naturally in terms of partial 
covering column changes. There is evidence for a relatively weak emission line from 
neutral iron with an equivalent width of only about 65 eV. After the flux transition, 
the strength of the edge feature is greatly reduced, suggesting a large reduction in the 
amount of partial covering. For a large region of statistically acceptable chi-squared 
parameter space, the luminosity of Cir X-l, after correction for partial covering, need 
not change during the transition. We discuss models for the partial covering and 
suggest that X-ray scattering by electrons may be important. Aspects of the Cir X-l 
spectrum are very similar to those of Seyfert 2 galaxies with Compton-thin tori. 
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1 INTRODUCTION 

Circinus X-l (hereafter Cir X-l) is a peculiar and poorly 
understood X-ray binary with an approximately 16.6 day 
period in which the compact object is thought to be a neu- 
tron star (Tennant, Fabian & Shafer 1986). The neutron star 
magnetic field strength, companion star spectral type and 
binary orbital parameters are not well determined. The lack 
of X-ray pulsations and the presence of type 1 X-ray bursts 
suggest the neutron star magnetic field strength is less than 
~ 10 11 G, and it has been argued that Cir X-l may be an 
'atoll' type X-ray binary with a magnetic field strength of 
less than ~ 10 9 G (Oosterbroek et al. 1995 and references 
therein). However, the broad band X-ray spectrum of Cir 
X-l has been argued to resemble that of a high magnetic 
field neutron star rather than that of a low magnetic field 
neutron star (Maisack et al. 1995; compare with Sunyaev 
et al. 1991), and if Cir X-l is associated with G321.9 - 0.3 



(see Stewart et al. 1993 and references therein) it would per- 
haps be difficult to understand such a weak magnetic field in 
such a young neutron star (see figure 1 of Phinney & Kulka- 
rni 1994). Type 1 bursting activity from Cir X-l appears 
to be highly intermittent, and to our knowledge no bursts 
have been seen either before or after the EXOSAT series of 
observations. The optical counterpart to Cir X-l has been 
unambiguously identified by Moneti (1992), but reddening 
uncertainties and potential accretion disc contributions have 
prevented a precise determination of the spectral type of the 
companion star. The orbit has been argued to be highly ec- 
centric due to the violent variability seen near zero phase 
(probably associated with the periastron passage of the neu- 
tron star), although proper orbital parameters are lacking. 
Given the orbital period and claimed space velocity of Cir 
X-l, it is probable that its orbit was at least initially highly 
eccentric (see figure 10 of Brandt & Podsiadlowski 1995). 

The iron K spectral features of Cir X-l and their rele- 
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vance to understanding the X-ray properties of Cir X-l more 
generally are also a subject of uncertainty. Very strong and 
potentially broad iron Ka emission lines have been discussed 
(e.g. Tennant 1985; Miyamoto & Kitamoto 1985; Makino 
1993; Gottwald et al. 1995; Maisack et al. 1995), but it is 
known that complicated iron K spectral features can be dif- 
ficult to interpret using data with low spectral resolution 
(e.g. Tennant 1988a; Gottwald et al. 1995). In particular, in 
Cir X-l models where a large amount of partial or total cov- 
ering is present (e.g. Ikegami 1986; section 3.2 of Tennant 
1988b; Inoue 1989; section 5 of Oosterbroek et al. 1995), one 
might expect the presence of a strong iron K edge. In this 
paper we use the excellent spectral resolution and sensitivity 
of the Japanese/USA Advanced Satellite for Cosmology and 
Astrophysics (hereafter ASCA; Tanaka, Inoue & Holt 1994) 
to probe the iron K complex of Cir X-l near one of its zero 
phase transitions. 

Following section 3 of Stewart et al. (1993) we take the 
distance to Cir X-l to be greater than 6.7 kpc. Neither the 
optical extinction nor the interstellar X-ray column to Cir 
X-l are precisely known, but the X-ray column is thought to 
be about (1-2) xlO 22 cm" 2 (e.g. Predehl & Schmitt 1995). 
We note that significantly smaller X-ray column values and 
arguments based on them (e.g. section 3 of Stewart et al. 
1991) may have trouble explaining the strength of the dust 
scattering halo of Cir X-l as shown in figure 7 of Predehl & 
Schmitt (1995). 

In some of the spectral fitting below we shall implic- 
itly assume Morrison & McCammon (1983) abundances (e.g. 
when we fit cold absorption and partial covering models) . To 
our knowledge there is no evidence which suggests that these 
abundances are incorrect, and scalings to different abun- 
dances may be made from our results. 



2 OBSERVATIONS AND DATA REDUCTION 

Cir X-l was observed with ASCA during AO-2 on 1994 Au- 
gust 4-5 using both Solid-state Imaging Spectrometer CCD 
detectors (SISO and SIS1) and both Gas Imaging Spectrome- 
ter scintillation proportional counters (GIS2 and GIS3). The 
energy resolutions of the SIS and GIS detectors are f» 2 per 
cent and w 8 per cent at the iron K complex. The observa- 
tions started at 20:13:09 UT on 1994 August 4 and ended 
at 10:30:45 UT on 1994 August 5. They covered the time 
of zero phase as given by equation 1 of Glass (1994) for 
N = 392. 

The SIS detectors were operated in 1 CCD mode, and 
the most well-calibrated SIS chips were used (chip 1 for SISO 
and chip 3 for SIS1). Bright mode was used for high bit rate 
data and fast mode was used for medium bit rate data. We 
shall use only the bright mode data in this paper, and the 
telemetry limit is 256 count s _1 SIS -1 for these data. The 
telemetry is saturated during the latter part of the observa- 
tion, and we only analyze the SIS data for which the teleme- 
try is not saturated. Cir X-l was placed at the nominal 1 
CCD mode pointing position. 

The GIS were operated in PH mode, and the bit as- 
signment was changed so as to give 6 bits of X and Y posi- 
tional information and an extra 4 bits of timing information 
(10-6-6-5-0-4 in the notation of section 3.3 of Day et al. 
1995b). We thus obtain 64 x 64 pixel images (i.e. RAWXBINS 



and RAWYBINS each have the value 64; note however that the 
GIS field of view is circular) and 1024 channel energy spec- 
tra in the 1-10 keV band. The telemetry limit is 128 and 
16 count s" 1 CIS" 1 for the high and medium telemetry bit 
rates, respectively. The telemetry was completely saturated 
for the medium bit rate data, and was saturated for the lat- 
ter part of the observation even at high bit rate. We shall 
use only the high bit rate data in this paper. When GIS 
photons are lost due to telemetry saturation, they are lost 
in a way that leaves the spectral shape unaffected. Thus GIS 
data that suffer from telemetry saturation can be analyzed 
reliably after suitable corrections (see below). 

We have used the 'Revision 1' processed data from God- 
dard Space Flight Center (GSFC) for the analysis below 
(Day et al. 1995a), and data reduction was performed using 
FTOOLS and xselect (see Day et al. 1995b for a descrip- 
tion of these packages and their application to ASCA data 
analysis). In this paper we focus on the iron K spectral fea- 
tures observed from a very bright source, and we have chosen 
our basic data screening criteria after consultation with the 
GSFC ASCA help service. The basic screening criteria wc 
adopt in this paper are listed in Table 1. We have used the 
SISCLEAN and GISCLEAN software in FTOOLS as appropriate, 
and for the SIS detectors we only use event grades 0, 2, 3 
and 4. 

The ANG_DIST housekeeping parameter is unacceptably 
large for the first « 1500 s of the observation as expected, 
but then it remains within the criteria of Table 1. 



3 ANALYSIS 

We first extracted images for Cir X-l in each detector and 
examined them. At low energies the profile of Cir X-l has 
extent over that expected for a point source due to its dust 
scattering halo (see Predehl & Schmitt 1995 for a detailed 
study of the dust scattering halo using higher spatial res- 
olution data from the ROSAT PSPC). This dust halo will 
seriously hinder X-ray imaging searches for jet emission. At 
higher energies (above ~ 5 keV), the profile of Cir X-l is 
consistent with that of a point source (this is as expected 
since dust scattering halo intensities drop strongly with in- 
creasing energy; see section 3.6 of Predehl & Schmitt 1995). 
There are no other strong X-ray sources in the fields of view. 
However, the GIS detectors also show weak apparent emis- 
sion to the south of Cir X-l at their edges (this emission is 
not from the calibration sources). This 'shoulder- like' emis- 
sion is part of the point spread function and is common 
in observations of bright sources (K. Ebisawa, private com- 
munication). We are confident that this emission originates 
from Cir X-l and is not from another X-ray source because 
its variability tracks the variability from Cir X-l described 
below (with reduced absolute amplitude). 

We use circular source cells centered on Cir X-l with 
radii of « 4.5 arcmin and » 9 arcmin for SIS and GIS 
source count extraction, respectively. The choice of back- 
ground regions is not straightforward for the SIS detectors. 
The brightness of Cir X-l requires us to use small circular 
regions of radius w 1.2 arcmin near the northeastern (SISO) 
and southeastern (SIS1) chip corners. Even with these re- 
gions, there is still some contribution from Cir X-l to the 
background cells. For the GIS detectors we use circular back- 
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Table 1. Data screening criteria for the SIS and GIS detectors. Details on these screening criteria may be found in section 5.2 of Day 
et al. (1995b). 



SIS screening criteria 


GIS screening criteria 


Description 


SAA= 


SAA= 


Only use data taken outside the South Atlantic Anomaly (boolean variable) 


ELV> 5 


ELV> 5 


Earth elevation angle for target (degrees) 


CQR> 4 


CQR> 4 


A/f itiimn m r"i i"f~— n~FF vicriHit^r ({~tr>\l f*~^-\ 

1V11111111U.111 V_. LI L> LJll lltllUlLy 1 vJJL, V ) 


ANG_DIST> 0.0 


ANG_DIST> 0.0 


Angular distance of field of view from specified direction (degrees) 


ANG_DIST< 0.01 


ANG_DIST< 0.01 


Angular distance of field of view from specified direction (degrees) 


T_SAA< or T_SAA> 16 




Avoid high background data just after SAA passage (seconds) 


S0_PIXL1> 




Number of events detected by chip 1 above the event threshold (SIS0 only) 


S0_PIXL1< 400 




Number of events detected by chip 1 above the event threshold (SIS0 only) 


S1_PIXL3> 




Number of events detected by chip 3 above the event threshold (SIS1 only) 


S1_PIXL3< 400 




Number of events detected by chip 3 above the event threshold (SIS1 only) 




G2_L1> 0.0 


Only use data inside the PHA discriminator range (GIS2 only) 




G3_L1> 0.0 


Only use data inside the PHA discriminator range (GIS3 only) 



ground cells with radii of ~ 2.0 arcmin that are essentially 
free of emission from Cir X-l. We have found that different 
background regions have very little effect on the results of 
the analysis presented below. 



3.1 Variability near zero phase 

3.1.1 Count rate variability 

Source and background light curves were extracted for all in- 
struments. We have corrected the GIS light curves for dead 
time effects by multiplying each point by the (time vary- 
ing) factor 1/(1-G2_DEADT) for GIS2 and 1/(1-G3_DEADT) 
for GIS3. The G2_DEADT and G3.DEADT parameters are ob- 
tained from the ASCA househeeping file, and they include 
the effects of both instrumental and telemetry dead time. 
When we quote GIS count rates below we implicitly include 
the correction for detector dead time. We found good gen- 
eral agreement between the shapes of the source light curves. 
There was also good general agreement between the shapes 
of these light curves and the shapes of the G2_L1 and G3_L1 
light curves. There is no evidence for strong background 
flaring, even when the screening criteria shown in Table 1 
are neglected (variability in the SIS background light curves 
can be attributed to the small contamination from Cir X-l 
itself). 

In Figure 1 we show the screened GIS 2 light curve of 
Cir X-l in the 1-10 keV band. Cir X-l is in a relatively low 
count rate state (w 30 count s^ 1 ) for about 6000 s, and it 
then makes an abrupt transition to a state where the mean 
GIS count rate is a factor of ~ 10 times higher. Given the 
detector position of Cir X-l, it has a GIS count rate of about 
0.05 times that of the Crab during the low count rate state 
and a GIS count rate of about 0.55 times that of the Crab 
during the high count rate state. The high count rate state 
shows strong and rapid variability for its first m 6000 s, 
but Cir X-l then settles into a calmer mode of behaviour. 
The Eddington limit of Cir X-l for hydrogen-poor material 
is thought to correspond to a flux of about 2.4 x 10~ 8 erg 
cm~ 2 s _1 (see section 2 of Stewart et al. 1991 but be wary of 
partial covering effects as described below), and in the high 
count rate state the flux after correction for only Galactic 
absorption is comparable to the Eddington limit flux (see 
Table 2 for numerical flux values). 



The basic shape of the zero phase light curve has been 
seen to change greatly over the years (see the discussion in 
section 3 of Tennant 1988b and Whitlock & Tyler 1994). Our 
observations do not show evidence for basic shape changes 
from the GINGA all-sky monitor and pointed observations 
reported by Tsunemi et al. (1989) and Makino (1993), al- 
though the rise after zero phase may be to a somewhat lower 
level H. We have verified this via further examination of the 
unpublished GINGA all-sky monitor data which show that 
Cir X-l does not appear to significantly exceed 1.3 Crab in 
the 1-20 keV band during 1987-1992 (note that the '1.8' and 
'0.6' in the caption for figure 8 of Tsunemi et al. 1989 are in- 
correctly transposed; see figure 7 of Tsunemi et al. 1989 and 
H. Tsumeni, private communication). RXTE all-sky moni- 
tor observations from 1996 show flux dips before zero phase 
as well as flux rises up to a remarkable 2-3 Crab shortly af- 
ter (within about 0.2-0.6 days) zero phase (see Shirey et al. 
1996). It is hard to determine whether the data from the 
GINGA all-sky monitor or the RXTE all-sky monitor are 
more germane for comparison with our 1994 ASCA data, 
although we do not see evidence for a count rate rise as dra- 
matic as is now seen by RXTE. In any case, we note that 
the observed variability supports the discussion in section 3 
of Glass (1994) regarding ephemerides. 

The high count rate state variability at 28000-36000 s 
in Figure 1 is of concern since the two rises appear to have 
similar shapes. We have manually examined the ANG_DIST 
parameter during this time, and while small fluctuations of 
w 9 arcsec or less are present, they do not appear to be 
able to explain the variability. This variability is also seen 
in the G2_L1, G2.LDHIT, G3_L1 and G3.LDHIT housekeeping 
parameters, and the G2_H0, G2_H2, G3_H0 and G3_H2 house- 
keeping parameters appear to be reasonably behaved during 



* We also made ASCA observations of Cir X-l on 1994 August 
26. These observations are about 5 days past zero phase according 
to equation 1 of Glass (1994) for N = 393, and the observed mean 
GIS2 count rate is m 350 count s —1 with 10 per cent variability. 
Comparison of this fact with Figure 1, figure 1 of Tennant (1988a), 
figure 8 of Tsunemi et al. (1989) and RXTE all-sky monitor data 
suggests that the 'decay' of the flux from Cir X-l after zero phase 
may be progressively occuring less rapidly. Further observations 
to search for such trends are needed. 
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Figure 1. Screened ASCA GIS2 light curve of Cir X-l after correction for detector dead time (see the text for the correction method 
used). The light curve is for the energy band 1.0—10.0 keV (GIS channels 85—849). The abscissa is in seconds after 20:37:12 UT on 1994 
August 4, and the data point bin size is 16 seconds. The vertical dotted line shows the predicted time of zero phase from equation 1 of 
Glass (1994) with N = 392. We shall refer to the first 6000 s of data as the 'low count rate state' and the data after this time as the 
'high count rate state.' There is very rapid and large amplitude variability between 6000-12000 s, and this variability causes the rough 
appearance of the data in this time interval. 



this time. The cut-off rigidity is well above 4 GeV c 1 during 
this time. 



3.1.2 Spectral variability 

In Figure 2 we show GIS 2 hardness ratios as a function of 
time. We define HR1 to be the ratio of the count rate in GIS2 
channels 256-849 to that in channels 85-255 (these energy 
ranges approximately correspond to 3.0-10.0 keV and 1.0- 
3.0 keV). We define HR2 to be the ratio of the count rate 
in GIS2 channels 171-255 to that in channels 85-170 (these 
energy ranges approximately correspond to 2.0-3.0 keV and 
1.0-2.0 keV). We define HR3 to be the ratio of the count 
rate in GIS2 channels 426-849 to that in channels 256-425 
(these energy ranges approximately correspond to 5.0-10.0 
keV and 3.0-5.0 keV). Note that we have reduced the length 
of the abscissa of Figure 2 relative to that of Figure 1 to only 
show the period of strongest spectral variability. 

During the low count rate state (the first 6000 s of Fig- 
ure 2) Cir X-l shows strong HR1 and HR3 variability but 
much less HR2 variability. The HR1 and HR3 variability of- 
ten has a flare-like character. The HR2 data during the low 
count rate state can be fit acceptably with a constant value 



of 0.969 {xl = 0.97 for 181 de grees of freedom), although 
there do appear to be some small systematic residuals. In 
light of this result we have examined the 1.0-3.0 keV light 
curve, and during the low count rate state the variability in 
this band is much less than that at higher energies. While 
the data are not strictly consistent with a constant model 
(xt = 1-17 f° r 181 degrees of freedom), no systematic count 
rate fluctuations by more than 20 per cent from the mean 
are apparent. In addition, the flare-like behaviour seen at 
higher energy is not obviously visible below 3 keV. 

All three hardness ratios are seen to change significantly 
after the transition from the low count rate state to the 
high count rate state, although the change at high energies 
appears to be less than that at low energies. 

After 12000 s (not shown in Figure 2), much less hard- 
ness ratio variability is seen. In particular, there is no ap- 
parent spectral variability corresponding to the count rate 
variability seen in Figure 1 between 28000-36000 s. Statis- 
tical analysis shows that there is a slow general decline of 
at least HR1 and HR3 with time. HR1 declines from 1.33 to 
1.22, and HR3 declines from 0.46 to 0.44. 
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Figure 2. ASCA GIS2 hardness ratios of Cir X-l as a function of time. The ordinate labels denote energy ranges in keV. The abscissa 
offset is the same as that for Figure 1, and we have only shown the start of the observation where the strongest spectral variability is 
observed. The data point bin size is 16 seconds. Note that HR3 can vary strongly while HR2 remains constant. 



3.2 Spectra near zero phase 

Our ASCA observation provides the highest spectral reso- 
lution view of Cir X-l yet possible and allows a sensitive 
probe of its iron K spectral features. To model the contin- 
uum underlying these features, we have used a two black- 
body model. As we shall show below, this model fits the con- 
tinuum in the ASCA band well even if the fitted blackbody 
parameters may not have precise translations into physical 
quantities. This continuum model was also found to fit the 
broader band («2-25 keV) and very high statistical quality 
GINGA data the best (Makino 1993). We have verified that 
our conclusions regarding the iron K spectral features are 
not materially altered if we use different continuum models 
such as Comptonized blackbodies (this is due to the good 
spectral resolution and good statistics). In addition, if we 
ignore the data below 3 keV and just use a power law to fit 
the 4-10 keV data our results are not qualitatively changed. 
In all fits below we include the effects of interstellar absorp- 
tion. In our fits with partial covering we, of course, take 
the interstellar absorption to be unaffected by the partial 
covering. 

Photon pile-up effects are a concern for ASCA SIS ob- 
servations of bright sources (see Ebisawa et al. 1996 for a de- 
tailed discussion of this issue). Cir X-l has a mean count rate 
during the low count rate state of w 30 count s -1 SIS -1 , and 



photon pile-up effects are unlikely to compromise the anal- 
ysis below. We have used ring-shaped extraction regions to 
verify that our results are not strongly influenced by photon 
pile-up effects. We note that the pile-up spectrum is smooth 
(see figure 11 of Ebisawa et al. 1996) and cannot falsely im- 
print the sharp spectral features we describe below. 

We use the SIS redistribution matrix files (rmf) from 
1994 Nov 9 and the CIS rmf from 1995 March 6. We generate 
our ancillary response files (arf) using the ASCAARF software. 
We have used the deadtime program in FTOOLS to correct 
for dead time effects. We model the X-ray spectra of Cir 
X-l using the X-ray spectral models in the xspec spectral 
fitting package (Shafer et al. 1991). 



3.2.1 Low count rate state spectral fitting 
Average spectral properties 

In order to obtain a 'first look' at the spectrum in the low 
count rate state, we have binned spectra for all four detectors 
using all of the low count rate state data. We group these 
spectra so that there are at least 20 photons per spectral 
data point (we shall implicitly adopt this grouping for all 
spectra below unless stated otherwise) . Figure 2 shows that 
there is significant spectral variability within the low count 
rate state, and thus these spectra must be regarded to be 
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Figure 3. Average SISO (solid circles) and SIS1 (open circles) spectra for the low count rate state data. The best fitting absorbed two 
blackbody model is shown along with the data-to-modcl ratio. Note the sharp drop in the data-to-model ratio just above 7 keV. The 
SIS detector resolution at the iron complex is about 140 eV, and the shown data points in this energy range are essentially statistically 
independent. The sharp drop can be interpreted as an iron edge from partial covering. 



crude and only meaningful in a time-averaged sense. How- 
ever, the excellent statistics afforded by this first look give 
useful insight. We note that an unlikely conspiracy would 
be required for spectral variability effects to produce sharp 
spectral features at the energies of known strong atomic fea- 
tures. 

In Figure 3 we show average SISO and SIS1 spectra for 
the entire low count rate state. The two spectra agree well 
in shape, and the GIS2 and GIS3 spectra are consistent with 
those in Figure 3 (they are not displayed only for reasons 
of clarity). We have used the two blackbody model with 
simple cold absorption (see above) to fit these data, and 
we show the ratio of the data to the best fitting model. In 
this fitting the blackbody normalizations are all taken to 
be free parameters, but the other spectral parameters are 
tied together (we adopt this fitting prescription for all fits 
below unless stated otherwise). The most noticeable feature 
in the ratio plot is the large and sharp drop just above 7 
keV. A natural interpretation for this drop is that it is due 
to an absorption edge from a large column of neutral or 
nearly- neutral iron atoms. The broad 5.3-7.1 keV excess and 
smaller 3.8-5.3 keV deficit may be understood as statistical 
artefacts from the fit trying to compensate for the presence 
of the large edge (fits which include the edge reduce the 
amplitudes of these features to where they are insignificant, 



see below). A large edge from iron would be expected in 
models for Cir X-l where a large amount of obscuring matter 
lies between the X-ray generation region and Earth, and one 
would not be naturally expected in other models. The iron 
structure could perhaps also be interpreted as broad iron 
line emission from matter moving with relativistic speed, 
and we shall critically examine this possibility below. There 
is also some secondary residual structure between 6-7 keV, 
although this structure is not as strong as that just above 7 
keV. There are no strong residuals in the 1-3 keV band. 

'Calm' set of data 

In order to probe the spectral properties of the low count 
rate state in a more precise manner, we have carefully se- 
lected a 450 s interval of data where there is relatively little 
flux and spectral variability. This interval begins at 21:55:31 
UT. We have binned spectra for all four detectors and fit 
these spectra using the two blackbody model with simple 
cold absorption. We derive the results shown in column 1 
of Table 2. This model can be ruled out with greater than 
94 per cent confidence, and similar systematic residuals to 
those shown in Figure 3 are present. 

The depth of the putative edge seen in the residuals and 
in Figure 3 is suggestive of intrinsic absorption by a column 
of ~ 10 24 cm~ 2 even though we do not see significant ab- 
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Table 2. Spectral fitting results. 

'CR' stands for 'count rate', 'PC stands for 'partial covering' and 'BB' stands for blackbody. 

The fitting models used below are defined in Shafcr ct al. (1991). When wc fit blackbody models we are using the 'bbodyrad' model, 
and when we fit partial covering models we are using the 'pcfabs' model. The blackbody model normalizations (Ai and A2) are formally 
related to the radii of the emission regions, although as commented in the text the derived radii are likely to only be characteristic 
quantities. A 'bbodyrad' normalization is formally equal to the square of the emission region radius in kilometers divided by the square 
of the distance to Cir X-l in units of 10 kpc. 

The errors for all fits are quoted for 90 per cent confidence, conservatively taking all free parameters to be of interest other than absolute 
normalization (Lampton, Margon & Bowyer 1976). 

Fluxes and normalizations are quoted for SISO when possible. When SIS data are not available due to telemetry saturation we quote fluxes 
and normalizations for GIS2. All fluxes and luminosities are for the 1—10 keV band. When we compute unabsorbed fluxes we remove the 
effects of both interstellar absorption and intrinsic partial covering. 

'd.o.f.' stands for 'degrees of freedom'. The P(x 2 I v ) chi-square probability function is defined in section 6.2 of Press et al. (1989). 



Fit 

parameter 
name 


'Calm' low 
CR state 
No PC 
(1) 


'Calm' low 
CR state 
PC on both BB 
(2) 


'Calm' low 
CR state 
PC on hot BB 
(3) 


HR3> 0.9 low 

CR state 
PC on both BB 
(4) 


High 
CR state 
No PC 

(5) 


High 
CR state 
PC on both BB 
(6) 


Interstellar 
N H /(10 22 cm" 2 ) 


1 77+0.11 
L -' ' -0.11 


1 vq+0-15 
1 - ' y -0.12 




1 09+0.23 
1 - 82 -0.20 


l-«2_ og 


1 70+0.11 


fcTi (keV) 


0.48±H 3 3 




o.48tg:E 


46+ 05 


0.59±H 3 3 


59+ 07 
U - OM -0.05 


Ai 


3710l 7 ?° 


(6.5±t;°) x io 4 


4150+ 1760 


(4.o± 2 ;°) x 10-4 


8950+ 2 ° 3 2 ° 


(1.5±l° 7 5 ) x 10 4 


kT 2 (keV) 


1.49i ;^ 


1 16+ - 19 

1 -' LO -0.15 


1 15+ - 17 
J - LO -0.15 


1 - ol -0.13 




1 QO+0.04 

L - 6 °-0.05 


A 2 


19l« 


1 10 -460 


o-in+2040 


201+ 120 


565tg 


1080^730° 


Partial covering 
percentage 




qq a+4-8 


qq s+ 4 - 7 
yd -°-10.8 






42 +54 

^ z -31 


Partial covering 
N H /(10 22 cm- 2 ) 




162i 4 4 ? 








19 o+109 


Absorption uncorrected 
flux/(10 — 9 erg cm - 2 s^ 1 ) 


1.5 


1.5 


1.5 


3.5 


24 


24 


Absorption corrected 
flux/(10- 9 erg cm- 2 s" 1 ) 


2.5 


38 


16 


25 


33 


41 


x 2 

d.o.f. 


1001 
932 
1.074 


931.6 
930 
1.001 


931.3 
930 
1.001 


1075 
1079 
0.996 


1389 
1279 
1.086 


1364 
1277 
1.068 


P(x 2 v) rejection 
probability (if significant) 


94.2 








98.3 


95.5 



sorption over the expected Galactic value at lower X-ray 
energies. This large differential column could either be ex- 
plained by partial covering from thick matter or by absorp- 
tion from ionized matter. Due to the fact that the putative 
edge energy corresponds to that of neutral or nearly-neutral 
iron, partial covering seems more plausible than warm ab- 
sorption. If we include partial covering of both blackbodies, 
we derive the results shown in column 2 of Table 2. The fit is 
statistically very good. In Figure 4 we show partial covering 
chi-squared contours for this fit as the solid set of contours. 
Note that a very large partial covering fraction is required 
by these data and that as a result the partial covering cor- 
rected isotropic luminosity for the low count rate state is 



of order the Eddington luminosity (it is difficult to be more 
specific than this because the observed X-ray bursts could 
have been affected at some level by partial covering too) . 

Narrow emission lines from iron have been seen in X-ray 
binaries (e.g. Ebisawa et al. 1996), and Figure 3 suggests the 
presence of an emission line from neutral or nearly-neutral 
iron. We have thus added a 6.40 keV (see House 1969) Gaus- 
sian emission line with a — 10 eV to our fit. This addition 
gives x 2 = 924.3 for 929 degrees of freedom, a significant 
improvement (the other fit parameters are not changed sig- 
nificantly by the addition of the line). Adding narrow lines 
at other energies (e.g. 4.0 keV, 6.1 keV, 6.6 keV) does not 
improve the fit significantly. The line flux is (8.9+5 g) x 10~ 3 
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Figure 4. Contours of the partial covering model for the 'calm' low count rate state data (solid curves near upper right), the low count 
rate state data with HR3 greater than 0.9 (dotted curves near upper left) and the high count rate state data (dashed curves). The partial 
covering is taken to affect both blackbodies (see the text). Contour levels for each set are for A\ 2 = 2.30 (68.3 per cent confidence for two 
parameters of interest), A\ 2 = 4.61 (90.0 per cent confidence for two parameters of interest) and A\ 2 = 9.24 (90.0 per cent confidence 
for five parameters of interest). Crosses show the best fit values. The low count rate state contours have been made using data from all 
four ASCA detectors, while the high count rate state countours have been made only using data from the CIS detectors. The partial 
covering luminosity correction increases towards the upper part of this diagram, especially for partial covering fractions near unity. 



photons cm~ 2 s _1 , and the line equivalent width is 65~^_4i eV 
(these errors are for A\ 2 = 2.71). The equivalent width we 
obtain from fitting the data shown in Figure 3 is compara- 
ble to this value. While there may also be small 6.7-6.9 keV 
lines from ionized iron, we are not able to formally prove 
the existence of such lines. The other detectors do not have 
?s 6.9 keV residuals as strong as the SIS0 « 6.9 keV point 
shown in Figure 3. 

We have also considered the possibility that only the 
higher temperature blackbody component suffers from par- 
tial covering. While figure 4 of Inoue (1989) appears to be 
at odds with this possibility, we note that HR2 does not 
vary during the low count rate state while HR3 does. The 
absence of HR2 variability, the absence of large 1.0-3.0 keV 
count rate variability (see above) and the presence of HR3 
variability could be suggestive of partial covering interposed 
between the higher and lower temperature blackbody emis- 
sion regions (but see our further discussion of this issue be- 
low in the context of partial covering column changes). From 
our fitting we derive the results shown in column 3 of Ta- 
ble 2. This is also a statistically good fit and thus from 
model fitting alone we cannot clearly determine whether 
the partial covering affects both blackbodies or only the 



higher temperature one. If we make partial covering con- 
tours such as those shown in Figure 4, they look very sim- 
ilar to the solid ones already shown. If we use the low and 
high temperature blackbody normalizations and a distance 
of 6.7 kpc to derive crude blackbody radii, the higher tem- 
perature blackbody has a radius of ^ 36 km and the lower 
temperature blackbody has a radius of ^ 52 km. Thus in 
this scenario the matter doing the partial covering would be 
roughly ~ 40 km from the neutron star. A very large den- 
sity of £ 10 22 cm~ 3 would be required for iron to be nearly 
neutral, as implied by the energy of the iron edge (see Ross 
1978 and Kallman & McCray 1982). This density is remark- 
ably high, being ~ 1400 times larger than the Eddington 
number density for accretion onto a 1.4 Mq neutron star 
(n E dd = I/Edd/47rr 2 c 3 m p = 1O 19 (M/M ) _1 cm~ 3 where 
r g = GM/c 2 ; see appendix C of Phinney 1983 and com- 
pare this density with section 5.6 of Frank, King & Rainc 
1992). The presence of such a dense partial coverer would 
be difficult to understand. Thus partial covering of just one 
blackbody, while statistically acceptable, does not appear to 
be physically plausible. The lack of strong Doppler blurring 
of the edge may also suggest that the partial coverer is sig- 
nificantly further than ~ 40 km from the neutron star. If the 
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partial covering affects both blackbodies, we derive a radius 
for the lower temperature blackbody of ~ 170 km, and the 
partial coverer could be far outside this radius. Thus sce- 
narios in which the partial covering affects both blackbodies 
significantly alleviate the density contrast problem. 

Now we consider the possibility that the iron K residuals 
are primarily due to a broad iron emission line from matter 
moving with relativistic speed near the neutron star (e.g. 
Brandt & Matt 1994 calculate the expected iron Ka line 
properties from matter in an ionized accretion disc around 
a neutron star). The blue horns of such lines can sometimes 
drop quite sharply at their high energy ends, and this could 
perhaps produce the sharp drop shown in Figure 3. We 
have fit the xspec 'diskline' model to the iron K spectral 
features. This model includes Doppler, transverse Doppler 
and gravitational redshift effects but does not include light 
bending (Fabian et al. 1989 describe the physics used in the 
'diskline' model but note that not all of the physics dis- 
cussed in Fabian et al. 1989 is actually implemented in the 
model). We constrain the line rest energy to be in the range 
6.40-6.97 keV (this is the iron K line energy range from 
House 1969 and Ross 1978). We fix the disc inner radius 
to be 6r g and constrain the disc outer radius to be larger 
than 7r g . We have found that the statistically best fits are 
obtained when we use the special '10' disc emissivity law 
appropriate for a black hole disc (see Shafer et al. 1991 and 
the XSPEC 'diskline' code). In this fit the line energy rises to 
6.97 keV and the inclination drops to degrees. The best fit- 
ting disc outer radius is 110r g , but this parameter is poorly 
constrained by these data. The line equivalent width is 440 
eV. This equivalent width, while large, is perhaps not im- 
possible for an ionized accretion disc (e.g. Matt, Fabian & 
Ross 1993). The fit gives \ 2 = 943.4 for 928 degrees of free- 
dom, and is thus statistically acceptable. From examination 
of the best fitting line model, it is apparent that the high 
energy drop of the blue horn of the line mimics the effect of 
an edge, as discussed above. We note that while a relativis- 
tic line cannot be formally falsified with these data, partial 
covering provides a slightly smaller \ 2 with fewer free pa- 
rameters (Ax 2 = 11.8 with 2 fewer parameters). In partic- 
ular, the drop just above 7 keV is uniquely and naturally 
specified in the partial covering model by the edge energy 
of cold iron. In addition, aspects of the observed spectral 
variability and low-to-high count rate state transition can 
be best explained in terms of partial covering (see below). 

In light of the results of this section, it appears likely 
that some earlier and lower spectral resolution observations 
which modeled the iron K complex of Cir X-l with only a 
line feature were not entirely correct. 

Other sets of data 

We have also fit other sets of data to investigate the extent 
to which our discussion above is generally applicable to the 
low count rate state. The first 500 s of data in Figure 1 
show relatively little flux and spectral variability, and the 
spectral residuals from these data are consistent with those 
discussed above. If we fit a partial covering model to these 
data, we obtain a set of contours that is very similar to the 
solid set shown in Figure 4. The partial covering fraction is 
again required to be greater than 80 per cent with very high 
statistical significance. 

Unfortunately, there are no other contiguous sets of low 



count rate state data during which both HR1 and HR3 re- 
main well-behaved for 200 s or more (200 s is about the 
time required to obtain a high quality spectrum during the 
low count rate state). In order to proceed further, we are 
forced to mollify our rigorous prescription of only looking at 
contiguous data sets with minimal spectral variability. We 
have selected a set of low count rate state data for analysis 
during which HR3 is larger than 0.9. The data used were 
chosen from that taken between 800-1700 s in Figure 2, and 
the total integration time was 280 s. Our choice of these data 
allows us to examine a significantly different HR3 range from 
that examined so far. In Figure 5 we illustrate the difference 
between a spectrum extracted from these data and one from 
the 'calm' data described above. As one would expect from 
Figure 2 and its accompanying discussion, the spectra agree 
well below 3 keV but are very different at higher energies. 
If, as argued above, the partial covering affects both black- 
bodies, the excellent agreement below 3 keV rules out any 
significant changes in the partial-covering fraction between 
the two data sets. However, the agreement does not rule out 
changes in the partial-covering column density between the 
two data sets as long as the column remains so high that it 
always entirely blocks the < 3 keV X-rays that strike it. In 
fact, the qualitative behaviour shown in Figure 5 is exactly 
what one would expect from a thick partial covering column 
that changes with time (this can be understood by exam- 
ining XSPEC partial covering models with partial covering 
columns between 5 x 10 23 cm" 2 and 10 25 cm" 2 ). If we fit 
the spectra from all four detectors with only a two black- 
body model, the fit is very poor and can be rejected with 
greater than 99 per cent confidence. If we add partial cov- 
ering, we obtain a good fit as shown in column 4 of Table 2. 
We note that the low energy spectral parameters (Nu, kTi 
and Ai) and the partial covering fraction agree well with 
what was derived for the 'calm' set of low count rate state 
data. While the kTi and A2 contours do not quite overlap 
with those derived for the 'calm' data, they are not strongly 
different and small discrepancies may be explained by the 
fact that we have fit a set of data which includes HR3 spec- 
tral variability. However, the most noticable change from the 
'calm' data set is the strong change in fitted partial cover- 
ing column. This can be clearly seen by comparison of the 
dotted and solid contours in Figure 4. The dotted contours 
are smaller than the solid contours because lower energy 
photons are now able to penetrate the column and thereby 
better constrain the fit. 

We have also examined other sets of data with HR3 be- 
tween 0.6-0.9. When we make partial covering contours for 
these data sets, the contours lie between the dotted and solid 
sets of contours in Figure 4. These data sets all have partial 
covering fractions larger than 0.8 with high statistical signif- 
icance. When we plot them in Figure 5 they lie in between 
the two curves shown (i.e. the spectrum 'unzips' from high 
to low energy as the partial covering column decreases). 

Physically, a large change in the partial covering column 
without some change in the partial covering fraction would 
appear somewhat surprising. Below we shall discuss ways in 
which the apparent partial covering column might change 
without an accompanying change in the apparent partial 
covering fraction. We shall also discuss how one might ob- 
serve apparent partial covering effects from an object as 
small as a neutron star and its inner accretion disc. 
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v Calm' and high HR3 spectra 
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Figure 5. SISO spectra extracted from the 'calm' set of data (lower spectrum) and a set of data which has HR3 greater than 0.9 (upper 
spectrum). Note that the spectra agree very well up to about 3 keV but then differ strongly at higher energies. The data sets have been 
fitted with double blackbody models which include partial covering (see the text). A partial-covering column change can explain the 
difference between these two spectra (see the text). 



3.2.2 High count rate state spectral fitting 

Due to telemetry saturation during the high count rate state, 
we have only been able to perform reliable analyses using 
the GIS detectors. To begin, we have extracted GIS spec- 
tra using only data taken within the 20000-24000 s inter- 
val of Figure 1. We group these spectra so that there are at 
least 40 photons per spectral data point. During this interval 
there is no strong spectral variability and the flux variabil- 
ity is rather small as well. If we fit the two blackbody model 
with simple cold absorption to these data, we obtain the 
results shown in column 5 of Table 2. While this fit is statis- 
tically unacceptable and can be rejected with greater than 
98 per cent confidence, examination of the data-to-model 
ratio reveals that fractional deviations from the model are 
much smaller in the high count rate state than the low count 
rate state. This important result is illustrated in Figure 6 
where we compare data-to-model ratios for the high and low 
count rate states. The difference in statistical fit quality is 
thus primarily due to the much larger number of photons 
available for spectrum construction during the high count 
rate state. Scrutiny of the high count rate state residuals re- 
quires caution due to the fact that any systematic residuals 
are so small that calibration uncertainties may have non- 
negligible effects on their interpretation (see appendices A 



and B of Ebisawa et al. 1996). Careful examination sug- 
gests that there may be a weak edge-like residual present 
around 7 keV, and this seems at least plausible in light of 
our results from the previous section. If we add partial cov- 
ering to our model as per the previous section, we obtain 
the results shown in column 6 of Table 2. The addition of 
partial covering gives a significant reduction of x 2 , although 
the fit is still not formally acceptable. There are no sys- 
tematic trends in the residuals that are obviously suggestive 
of spectral features, and we suspect that calibration uncer- 
tainties are contributing significantly to the x 2 value due to 
the very good statistics. If we include 2 per cent calibration 
uncertainties when we bin our spectra (see Tashiro et al. 
1995 for a discussion of GIS calibration uncertainties), we 
are able to obtain statistically acceptable x values. In Fig- 
ure 4 we show partial covering confidence contours as the 
dashed set of contours. The high count rate state contours 
in this diagram are larger due to the fact that the strength 
of the residuals that constrain the partial covering are much 
smaller (see Figure 6). We note that GIS calibration un- 
certainties almost certainly cannot make the low and high 
count rate state contours in Figure 4 consistent with each 
other (as is to be expected from Figure 6) . 

We have examined whether a narrow 6.4 keV emission 
line from iron is present in the high count rate state data 
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Figure 6. GIS2 data-to-model ratio plots for the low count rate state data (upper panel) and the high count rate state data (lower 
panel). Note that the abscissae and ordinates of both panels have the same scales, and that the low count rate state data have much 
stronger iron K complex residual structure than the high count rate state data. The spectral resolution of the data in this figure is about 
4 times less than that in Figure 3. 



by adding a Gaussian emission line model. The addition of 
the line is not statistically significant, and the best fitting 
line normalization is zero. The line normalization can be 
constrained to be less than 4.5 x 1CT 3 photons cm~ 2 s _1 
with 90 per cent confidence (A^ 2 = 2.71). This corresponds 
to a line equivalent width of less than 10 eV. While we find 
that the line equivalent width is much smaller than for the 
low count rate state, we cannot formally prove that the line 
normalization is much smaller. 

We have analyzed other sets of data taken after 12000 
s in Figure 1, and the iron edge appears to be generally 
much weaker during this time than during the low count 
rate state. In particular, spectra created from data taken 
during the two rises between 28000-36000 s in Figure 1 do 
not show iron K residuals any stronger than in the bottom 
panel of Figure 6. 



3.2.3 Spectral ratios 

As discussed by Inoue (1989), the ratio of two spectra that 
differ only in the fraction of thick partial covering imposed 
on them gives useful information about this partial covering. 
At X-ray energies of < 3 keV, where photoelectric absorp- 
tion blocks all X-rays striking the partial covering material, 
the spectral ratio should approach a constant value. This 



constant value gives the relative partial covering fractions in 
the two spectra. At higher X-ray energies, as photoelectric 
absorption becomes progressively less important, the spec- 
tral ratio should approach unity (this effect will not be ob- 
servable in the ASCA band for columns significantly larger 
than ~ ov^ 1 ~ 1.5 x 10 24 cm -2 ). Exactly this behaviour was 
reported by Inoue (1989), and we have made ratios of our 
spectra to look for it. 

We have made a ratio of the GIS2 spectrum from the 
'calm' set of low count rate state data to the GIS2 spec- 
trum from the bottom panel of Figure 6 ('spectral ratio 1' 
or 'SRI') as well as a ratio of the GIS2 spectrum from the 
HR3> 0.9 set of low count rate state data to the GIS2 spec- 
trum from the bottom panel of Figure 6 ('spectral ratio 2' 
or 'SR2'). We show a plot of these ratios in Figure 7. SR2 
starts to approach unity above 4 keV as its numerator's col- 
umn is progressively penetrated by higher energy X-rays, 
and SRI also starts to rise at the high energy end of the 
ASCA band (this behaviour will be able to be probed much 
better with RXTE). The < 3 keV behaviours of SRI and 
SR2 are quite similar. However, these behaviours are both 
significantly different from that shown in Figure 4 of In- 
oue (1989). They resemble more closely the behaviour seen 
by Ikegami (1986). Neither spectral ratio is approximately 
constant below 3 keV as would be expected under the hy- 
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pothesis of only a partial covering fraction change between 
the low count rate state and the high count rate state. 

Using the two blackbody continuum model with partial 
covering, we have also determined that the combination of a 
partial covering fraction change and a partial covering col- 
umn change cannot, by itself, explain the difference between 
the low and high count rate state spectra. We have done this 
by jointly fitting the GIS spectra for the 'calm' low count 
rate state data and the high count rate state data. In this fit- 
ting we tie all model parameters together other than those of 
the partial covering. The GIS2 and GIS3 blackbody normal- 
izations are separately tied together for each pair of spectra. 
We include 2 per cent calibration uncertainties. The result- 
ing joint fit is statistically unacceptable and can be rejected 
with greater than 99 per cent confidence. As expected from 
Figure 7, examination of the residuals shows that the low 
count rate state spectra are softer below « 2 keV than the 
high count rate state spectra. If we add a third blackbody 
that is unaffected by partial covering to our joint fitting, we 
can obtain an acceptable joint fit. The third blackbody has 
a temperature of about 0.5 keV but is poorly constrained. 

The fits of the previous paragraph and their residuals 
could signify either (1) a true change of the blackbody pa- 
rameters between the low and high count rate states (as 
well as a partial covering change) or (2) the presence of 
an additional soft component that is unaffected by the par- 
tial covering so that it can be seen during the low count 
rate state but is diluted away from observability during the 
high count rate state. As is clear from previous sections, 
our spectral fits alone cannot formally falsify the first pos- 
sibility. However, given the observed spatial extent at low 
energy the additional soft component of the second possibil- 
ity can naturally be explained as time-delayed X-rays that 
have been dust scattered into the line of sight by the inter- 
stellar medium. For a source at the distance and position 
of Cir X-l, the time delay is ~ 10 days (see equation 20 of 
Mauche &: Gorenstein 1986), a significant fraction of the or- 
bital period. Replacing the third blackbody of the previous 
paragraph by a model consisting of the first two blackbodies 
multiplied by E~ 2 does give a fair fit to the data so indi- 
cating that dust scattering is a plausible hypothesis. Our 
relatively short span of low count rate state data does not 
allow us to probe sensitively for temporal changes of the soft 
component due to the decay of the dust halo. 



4 DISCUSSION AND SUMMARY 

4.1 Partial covering: a direct flux plus electron 
scattered flux model 

Our observations strongly argue that partial covering is a 
crucial determinant of the behaviour of Cir X-l. We clearly 
see the spectral features expected from partial covering (Fig- 
ure 3), and the observed continuum spectral variability is 
also best explained in terms of it (Figure 5). These two sep- 
arate effects in conjunction are compelling. As mentioned in 
Section 3, a large change in the partial covering column with- 
out a corresponding change in the partial covering fraction 
(see Figure 4) is physically somewhat surprising. However, 
as we describe below, such a change could easily occur if 
we observed both direct, but attenuated, X-rays as well as 



Spectral ratios 




Energy (keV) 

Figure 7. GIS2 spectral ratios as labeled in the diagram ('CR' 
stands for 'count rate' and 'SR' stands for spectral ratio). 

X-rays that are electron scattered around the attenuating 
matter (this electron scattered component is local to Cir X- 

1 and totally distinct from the dust scattered component of 
the previous section). Cir X-l would then be like a Seyfert 

2 galaxy with a Compton-thin and time- variable torus, and 
the fitted partial covering would really represent total cover- 
ing plus a significant electron scattered component. Partial 
covering models have been used to fit Seyfert 2 X-ray spectra 
(e.g. Iwasawa et al. 1994). As will become clear, the direct 
flux plus electron scattered flux scenario obviates the geo- 
metrically difficult requirement to partially cover an object 
as small as a neutron star and its inner disc. 

In this scenario, the direct X-rays would be completely 
blocked at low energies but would penetrate through the 
absorbing matter above ~3-5 keV and thereby imprint the 
observed iron edge on the spectrum. The absorbing material 
could be the outer bulge of an accretion disc if the disc is 
seen nearly edge-on, and the column changes could occur as 
the disc rotates (dynamic changes of the partial coverer are 
possible as well) . For a disc outer radius of ~ 1 solar radius 
the disc rotation timescale is such that this scenario appears 
plausible. Several Seyfert 2 galaxies with Compton-thin tori 
also show large edges from neutral iron (e.g. Awaki et al. 
1991). 

The electron scattering of X-rays could arise in an accre- 
tion disc corona. One potential prediction of models with a 
large amount of electron scattering is that the disc flux from 
Cir X-l may be polarized over the expected interstellar value 
at wavelengths where the direct emission is highly absorbed 
(of course, one will also have to take into account polariza- 
tion dilution by the companion at wavelengths where this is 
relevant). The polarization fraction should change with time 
as the absorbing column, and hence the direct contribution, 
changes. If scattered photons suffered multiple scatterings 
the polarization could be significantly reduced (indeed, if 
the Thomson depth of the scattering 'mirror' were greatly 
below unity the implied isotropic luminosity would be pro- 
hibitive). If the jets of Cir X-l are perpendicular to its disc, 
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one might expect the polarization position angle to be per- 
pendicular to the jet position angle. 

From the previous paragraphs and observations of 
Seyfert 2 galaxies, one might expect to observe very large 
equivalent width (~ 1 keV) iron K lines from Cir X-l during 
the low count rate state. However, it is important to remem- 
ber that some direct continuum X-rays are able to penetrate 
through the obscuring matter at the energies of the iron K 
lines, and these will reduce the line equivalent widths. Fur- 
thermore, the solid angle subtended by the accretion disc 
bulge need not be large, and this would reduce the equiva- 
lent width of the iron line from neutral matter. Performing a 
calculation similar to that shown in figure 1.11 of White, Na- 
gase & Parmar (1995), we find that A£l/2n must be about 
0.1 to explain the iron line equivalent width, which is con- 
sistent with what one might expect from the outer bulge of 
an accretion disc. The ionized iron lines need not be strong 
if iron atoms in the scattering mirror are fully stripped of 
their electrons, and the large luminosity of Cir X-l allows 
iron atom stripping out to radii of ~ 0.2Rq or more (we 
take the Thomson depth to be r w 1 and require the ion- 
ization parameter to be £ ^ 10 4 erg cm s _1 ; see Kallman & 
McCray 1982). Lower atomic number atoms would also be 
fully stripped and thus would not produce lines. 

Finally, one might expect to see a broad iron line and/or 
smeared edge due to X-ray reflection off the surface of the 
inner accretion disc (e.g. Ross, Fabian & Brandt 1996 and 
references therein). No such features are clearly detected 
(see the bottom panel of Figure 6). There are three possible 
reasons for the absence of such features: (1) the magneto- 
sphere or an instability disrupts the inner accretion disc so 
little disc reflection takes place, (2) the disc inclination is 
high so that the reflection component is relatively weak (see 
George & Fabian 1991; note this is what we have suggested 
above) and (3) electron scattering in the corona, if its tem- 
perature is ^ 1 keV, renders features indistinguishable by 
Doppler broadening (such Doppler broadening has only a 
small effect on the continuum shape) . 

4.2 Partial covering and the abrupt count rate 
change 

Our spectral fitting above suggests a strong reduction in the 
amount of partial covering during the low to high count rate 
state transition (see Figure 1 and Figure 6). Furthermore, 
for our best fits, the partial covering corrected isotropic lu- 
minosities (hereafter 'the true luminosities') in the low and 
high count rate states are quite similar (the absorption cor- 
rected fluxes are 3.8 x 10 -8 erg cm -2 s _1 and 4.1 x 10~ 8 erg 
cm~ s~ , respectively). Indeed, we find that for a large re- 
gion of statistically acceptable chi-squared parameter space 
the true luminosity of Cir X-l need not change during the 
observed transition. It is more difficult to prove the stricter 
statement that the true luminosity of Cir X-l does not 
change, due to the fact that for very thick partial cover- 
ing the true luminosity and partial covering percentage are 
highly covariant quantities (cf. the contours in Figure 4). If 
the true luminosity of Cir X-l does not change significantly 
between the low and high count rate states, Cir X-l would 
have been constantly accreting at or perhaps even above the 
Eddington rate throughout our observation, despite the fact 
that its count rate varied by about an order of magnitude. 



Furthermore, as discussed in section 3.2.3, the underlying 
spectral shape need not change either provided there is a 
third soft component that is diluted to invisibility once the 
partial covering is reduced (again it is hard to prove that the 
underlying spectral shape does not change). We note that, 
despite the lack of strong low energy absorption changes, 
much of the behaviour studied in Dower, Bradt and Morgan 
(1982) could also result from partial covering changes pro- 
vided the partial coverer is so thick that essentially all the 
X-rays striking it are entirely absorbed. Indeed, the famous 
figure 4 of Dower et al. (1982) is highly suggestive of partial 
covering by very thick matter. 

However, as discussed by Dower et al. (1982) and other 
authors, the strong radio and infrared outbursts near zero 
phase cannot be obviously understood solely in terms of par- 
tial covering changes. Furthermore, due to the observed flux 
rises during the X-ray bursts, we know that the true lumi- 
nosity of Cir X-l cannot be the Eddington one at all times. 
The recent flux rises to 2-3 Crab seen in RXTE all-sky mon- 
itor data also suggest true luminosity changes. The detailed 
situation is likely to be very complex, although our obser- 
vations convincingly demonstrate that some form of partial 
covering is a key element. It is perhaps possible that the 
obscuring matter is ejected from the vicinity of the neutron 
star at the time of zero phase, perhaps to make the radio 
nebula or jets. Our iron line measurements are consistent 
with an ejection of matter although the statistics are not 
good enough to prove that such an ejection must occur (i.e. 
we cannot show that the line normalization must drop dra- 
matically). 

Our observations do not appear to give direct insight 
into why Cir X-l accretes at such a remarkable rate, and 
we have little to add to the theoretical speculation on this 
issue. It may well be difficult to gain further insight into this 
matter until the orbital parameters, companion star spectral 
type, proper motion and neutron star magnetic field of Cir 
X-l are better determined. Observations to measure these 
crucial system parameters are needed. 
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